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Observations made during the last ten years with the Chandra X-ray 
Observatory have shed much light on the cooling gas in the cen- 
ters of clusters of galaxies and the role of active galactic nucleus 
(AGN) heating. Cooling of the hot intracluster medium in cluster 
centers can feed the supermassive black holes found in the nuclei of 
the dominant cluster galaxies leading to AGN outbursts which can 
reheat the gas, suppressing cooling and large amounts of star forma- 
tion. AGN heating can come in the form of shocks, buoyantly rising 
bubbles that have been inflated by radio lobes, and the dissipation 
of sound waves. 

clusters of galaxies | active galactic nuclei | X-ray observations | Chandra X- 
ray Observatory 

Introduction 

Chandra observations of cavities in cluster centers in the 
last decade have brought the issue of feedback to the cen- 
ter stage in astrophysics. Feedback from AGN driven by su- 
permassive black holes in the cores of massive galaxies is now 
seen as a necessary ingredient to adequately describe the for- 
mation and evolution of galaxies and the clusters in which they 
sometimes reside. The energy injected by AGN is required to 
produce a galaxy luminosity function in simulations (e.g. [T]) 
that matches observations. This feedback may set the upper 
limit to the observed masses of galaxies and contributes to 
cluster preheating and the observed "entropy floors" seen in 
cluster centers. It may affect scaling relations, such as the 
relationship between X-ray luminosity and temperature and 
potentially affects cluster properties that are used for con- 
straining cosmological models, such as the gas mass fraction 
(e.g. [5]). Data from the Chandra X-ray Observatory over the 
last ten years reveal AGN at work in the centers of clusters, 
inflating bubbles that rise buoyantly through the intraclus- 
ter medium (ICM), sometimes producing shocks and sound 
waves. 

Many of the most spectacular Chandra images of AGN 
feedback are of interactions in the centers of cooling flow 
(sometimes referred to as "cool core") clusters of galaxies. 
Under the assumption of hydrostatic equilibrium, the X-ray- 
emitting gas in a cluster is the most dense at the cluster cen- 
ter. The gas will cool as it radiates, and the cooling time 
is inversely proportional to the density of the gas. There- 
fore, gas will cool at cluster centers before it cools at larger 
radii. If the cooling time is shorter than the cluster age, the 
central gas will cool, lose pressure, and outer gas must flow 
in to maintain hydrostatic equilibrium (see [3| for a review). 
Data from earlier X-ray observatories such as Einstein [4], the 
Roentgen SATellite (ROSAT) [5[, and the Advanced Satel- 
lite for Cosmology and Astrophysics (ASCA) jo, were used to 
give estimates of the gas inflow rates of typically hundreds 
or thousands of solar masses per year. The classic "cooling 
flow problem" has been that sufficient quantities of cool gas 
and/or star formation were not detected to match the predic- 
tions from the X-ray data. 

Chandra and XMM-Newton observations have shown that 
while the gas temperature drops in the centers of cooling 



flow clusters, the temperature reaches a threshold value below 
which there is little cooling. This is seen in the temperature 
proflles as well as high resolution spectroscopy. An important 
early result from XMM-Newton high resolution spectra from 
cooling flows was that the emission lines from cool gas were 
not present at the expected levels, and the spectra were well- 
fltted by a cooling flow model with a low temperature cutoff 
[T] [8]. This cutoff is typically one-half to one-third of the av- 
erage cluster temperature, and is usually not much below 1 
keV. The gas is either then cooling non-radiatively or being 
heated to higher temperatures. Throughout this paper, we 
use the term "cooling flow" to indicate clusters with central 
drops in temperature and large central density peaks above 
that of a typical /3-model. 

There is a clear association with the presence of a cooling 
flow in a cluster and the central galaxy hosting a radio source. 
Using Einstein data, [5] found « 70% of cooling flow clusters 
had giant elliptical galaxies that emitted in the radio, as com- 
pared to only « 20% of central galaxies in non-cooling flows. 
Recent studies, depending on the sample, find a fraction of up 
to 100% for cooling fiows and 45% for non-cooling flows [TD] 
using Chandra data. This association is consistent with the 
cooling gas supplying fuel to the central black hole, leading 
to the AGN outburst observed in the radio. As seen in Chan- 
dra images, the radio jets and lobes associated with the AGN 
signiflcantly affect the X-ray gas, inflating "cavities" or "bub- 
bles" in the cluster centers. Clear correspondence between 
regions of radio emission and deflcits in the X-ray are seen in 
a number of cases (e.g. Hydra A [111 I12j . Perseus [13], Abell 
2052, Abefl 2597 [l5l[16], Abefl 2199 [l^, Centaurus [18], 
AbeU 4059 JcF, '201, Abefl 262 [211 [22], MS0735. 6-^7321 [23], 
M87/Virgo 24 ). The X-ray gas, in turn, affects the radio 
lobes by confining and distorting them. Energy can be trans- 
ported to the ICM through bubbles that rise buoyantly to 
larger cluster radii, as well as through weak shocks and sound 
waves generated by AGN outbursts. The energy transferred 
to the cluster gas can diminish the cooling, giving an expla- 
nation for the lack of copious quantities of cool gas and star 
formation seen in the centers of cooling flow clusters. Com- 
posite images of A2052 and A262 are shown in Figs. 1 and 
2, respectively, where blue is the X-ray emission as seen with 
Chandra, red is optical emission, and yellow/green (in the 
cluster centers) is 1.4 GHz radio emission. 
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Cavity properties 

The sizes of the cavities or bubbles created by the radio lobes 
evacuating regions of the ICM vary widely from a few kpc (e.g. 
Abell 262 [HI [22]) to hundreds of kpc (e.g. MS0735. 6+7321 
[23]) in diameter. The cavities are often surrounded by bright 
rims or shells, where the gas from the cavity interiors has been 
compressed. In addition to the visual impression that the ra- 
dio lobes are removing the X-ray gas from the cavities, we can 
quantitatively confirm this by measuring the mass of gas in 
the shells and comparing it to the mass of gas that should fill 
the central region of the cluster estimated by extrapolating 
the density profile of the gas at outer radii inward. In Abell 
2052 |14l I25| , for instance, the mass in the shells compressed 
by the two radio lobes is approximately 5 x 10^" M©, consis- 
tent with the mass predicted to have filled the cavities, based 
on the extrapolation of the density profile. 

The X-ray gas in the compressed shells surrounding the 
radio lobes is generally found to be cool, and does not show 
evidence of strong shock heating as had been predicted with 
early models of radio source / ICM interaction (e.g. [26j). 
Cool gas is seen surrounding the inner radio lobes of, for ex- 
ample. Hydra A [27], Perseus A [28], 3C 317 in Abell 2052 
[H [25], and B2 0149-^35 in Abell 262 [21]. The cooling time 
within the shells is typically longer than the radio source age, 
suggesting that the gas in the shells had to cool closer in to the 
cluster center before being pushed out by the radio lobes. In 
the case of Abell 2052, the shell cooling time is approximately 
3 X 10* yr while the age of the radio source is approximately 
10^ yr. The bright X-ray shells are sometimes coincident with 
optical line emission, such as Hot, representing gas that has 
cooled to temperatures of approximately 10* K [141 1211 129| . 
Therefore, while most of the cluster gas cools to only one- 
third to one-half of the average ICM temperature, at least 
some quantity of gas can cool to much lower temperatures 
and sometimes results in star formation [30| . 

Assuming some geometry for the shells surrounding the 
radio lobes and fitting the spectra from these regions, the 
densities and pressures can be determined. Alternatively, the 
pressure in the X-ray gas is determined from undisturbed re- 
gions in the cluster at the same radii as the cavities, yielding 
similar results The pressures found are often higher (by 
approximately an order of magnitude) than the pressures de- 
rived from radio observations, assuming equipartition, for the 
cavities. The cavities would be expected to collapse quickly if 
they were indeed underpressured, and it is therefore assumed 
that some additional source of pressure must be filling the cav- 
ities. The origin of this additional pressure is still unknown, 
but candidates are low energy, relativistic electrons, a higher 
ratio of ions to electrons than is typically assumed, a magnetic 
field with a value higher than the equipartition value, or very 
hot, diffuse, thermal gas. Hot, diffuse, thermal gas has not 
been detected filling the radio cavities, but limits have been 
placed in a few cases (> 15 keV in Hydra A [27], > 11 keV in 
Perseus [28], and > 20 keV in Abell 2052 [25]). 

The cavities are eventually expected to break up due to 
Rayleigh- Taylor and Kelvin-Helmholtz instabilities. Rayleigh- 
Taylor instabilities are a possible explanation for the spur of 
gas jutting into the N cavity in Abell 2052 14, 25, 32_. The 
cavities can be maintained against these instabilities by mag- 
netic draping [33] or pressure related to bubble inflation [34| . 



Buoyant bubbles 

Since the density inside the radio cavities is much lower than 
that of the ambient gas, the cavities should be buoyant, and 
will rise out to larger cluster radii. Such outer cavities are 



called "ghost cavities" or "ghost bubbles" because they are 
thought to be the result of earlier AGN outbursts. In addition, 
if these cavities are filled with radio lobes, the lobes often have 
steep spectra, and radio emission at the commonly-observed 
radio frequency of 1.4 GHz is sometimes faint or absent. Radio 
observations at lower frequencies have shown that the radio 
sources extend into these outer cavities in many cases. 

Examples of systems with outer cavities are Perseus |13j . 
Abell 2597 [15] [16], Abell 262 [22], Abell 2052 [32], and Hydra 
A [12]. In all of these cases, when the cavities seemed devoid 
of radio emission at higher frequencies, lower frequency radio 
emission was found to extend farther from the cluster center 
to fill the outer cavities. 

Entralnment of cool gas. Cool gas can be entrained along with 
radio lobes as they propagate into the ICM. There is a clear 
association with arcs of cool gas and the radio jets and lobes in 
the M87/Virgo system [24U35l[36l . This is also seen in the cen- 
ter of Abell 133, where a cool X-ray plume appears related to 
a detached radio lobe or relic [37]. In the Perseus cluster [29], 
plumes of cool, ~ 10* K gas seen in Ha appear to have risen 
away from the cluster center, having been dragged by one of 
the outer cavities. The metallicity structure is also consistent 
with the entrainment scenario [38] . Analysis of deep Chandra 
observations of the Hydra A cluster reveals higher metallicity 
gas along the radio jets and lobes, indicating that the gas is 
being redistributed from the center (where the gas metallicity 
is typically higher) to larger radii in the cluster [38| . 

X-ray cavities as radio calorimeters. Rather than relying on 
radio observations to estimate the energy input from the ra- 
dio sources where the ratio of radiative power to kinetic power 
can vary widely, we can instead use observations of the effects 
that the radio sources have on the X-ray-emitting ICM to de- 
termine how much energy is being deposited into the cluster 
gas. As in [11] and [39], the energy input is a combination of 
the internal energy in the cavity and the work done in inflating 
the cavity, or the enthalpy: 

(7^^^ + ^^^-(T^^^ 

where P is the pressure measured in the X-ray shells (which 
are assumed to be in pressure equilibrium with the cavities), 
and V is the volume of the cavities. The value of 7 is the 
mean adiabatic index of the gas filling the cavities and is 4/3 
for relativistic gas or 5/3 for non-relativistic, monatomic gas. 
In the relativisitic case, the total energy input is then 4PV , 
and in the non-relativistic case it is (5/2)PV^. To estimate an 
average energy input rate, we need the radio source outburst 
rate. In some cases, this can be measured directly from the X- 
ray data, when there are multiple sets of bubbles. Taking the 
projected distances between outer cavities and inner cavities 
and calculating the time to rise buoyantly to the outer posi- 
tion, the time between outbursts can be estimated. The time 
between outbursts is typically a few tens of millions of years 
(e.g. Perseus [l3]; A262 [22, A2052 [32]). In cases without 
outer cavities, an approximate value based on other systems 
is assumed. 

Does energy input from bouyantly rising bubbles offset cool- 
ing?. We can compare the energy injection rate from AGN- 
inflated cavities that rise buoyantly in clusters to the lumi- 
nosity of cooling gas in the cluster centers. This cooling lumi- 
nosity is gvien by 

Lcooi = ^ — kT [2] 

2 flTTlp 
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where M is the mass deposition rate into the cluster center of 
the cooling gas, T is the upper temperature from which the 
gas is cooling, and ^ is the mean molecular weight with 
being the average mass per particle in the gas. In many cases, 
the energy input rate from the buoyant bubbles is enough to 
offset the cooling gas (e.g. AbeU 2052 [T11[2S1|32| with a rate 
of « 3 X 10*^ erg s"^ and Hydra A [TTl |27l O] with a rate 
of ~ 3 X 10''* erg s~'). In some cases, the bubbles did not, 
at least initially, appear to supply enough energy to the ICM 
(e.g. Abell 262 [2T], where the energy input rate is approxi- 
mately an order of magnitude too low). 

In addition to detailed studies of individual objects, stud- 
ies of samples of cooling flow clusters with central AGN have 
been undertaken gT] |^ gSl IH gS]. Thirty-three systems 
were included in a sample in g^, and it was found that, on 
average, APV of enthalpy per cavity supplies enough energy 
to substantially suppress cooling in clusters. 



Shocks 

Simulations of AGN outbursts in cluster centers show the cre- 
ation of shocks as well as cavities related to the radio lobes 
(e.g. gB]). These "cocoon" shocks can heat the cluster gas 
not only in the direction parallel to the propagation of the ra- 
dio lobes, but perpendicular to their propagation axis, as well. 
Therefore, energy can be distributed throughout the cluster 
center if shocks are present. 

Until recently, shocks were not seen associated with AGN 
radio lobes in cluster cooling flows, but in the last few years, 
several have been detected. A spectacular example is the clus- 
ter MS0735. 6-1- 7321 [23] which hosts extremely large cavities 
with diameters of approximately 200 kpc. The total energy 
injection required to inflate the cavities and produce the ob- 
served shocks is 6 X 10^' erg, making this the most powerful 
radio outburst known. The shock energy input rate is given 
by 

where P is the pre-shock pressure, 7 = 5/3, and 2-k /ui is the 
time interval between shocks [31]. A deep Chandra image of 
M87/ Virgo revealed an approximately spherical shock in the 
cluster center [24]. The shock is most clearly seen in a hard- 
band image, and is fairly weak with a Mach number of 1.2. 
A Mach 1.65 shock was found in the Hercules cluster [47| . 
with the total energy deposited 3 x 10^^ erg. There is a weak 
shock (Mach 1.3) outside of one of the bubble rims in the 
center of the Perseus cluster [48]. As is common, the density 
jump measured is consistent with a shock model, but a rise 
in temperature is not measured. Two shocks are seen sur- 
rounding the center of AbeU 2052 32;, both with Mach 1.2. 
The temperatures across the shocks are consistent, within the 
errors, to the jumps expected in the shock models, but the 
best-fitting temperature values are constant across the fronts. 
If an increase in temperature is not measured behind a shock, 
[49j has suggested that shocks may be isothermal if conduc- 
tion is efficient in these regions. However, the temperature 
rises are very difficult to detect given the narrow width of the 
shock fronts, adiabatic expansion, and the difficulty in mod- 
eling the large amount of cluster gas at larger radii projected 
onto these regions when fitting the X-ray spectra. Based on 
further analysis of the Perseus system, [48] conclude that mix- 
ing of post-shock gas with cool gas associated with optical-line 
emission may result in the lack of temperature rise seen asso- 
ciated with the shock. 



Sound waves 

In addition to weak shocks in the centers of cooling flow clus- 
ters, sound or pressure waves have also been observed. These 
features are most prominent in the Perseus cluster [491 150] . 
They appear as 5 — 10% increases in pressure in ripple-like 
patterns surrounding cluster centers. Changes in tempera- 
ture are not seen in these regions, and it is likely that their 
energy is dissipated by viscosity. There may be one or more 
set of ripples corresponding to each AGN outburst depending 
on variations in energy output that can occur during a single 
outburst. Assuming they are associated with single outbursts, 
the ripple separations in Perseus imply a repetition rate of 10^ 
yr [50]. As shocks propagate outward from cluster centers, 
they weaken and may appear as sound waves. In addition to 
the Perseus cluster, sound wave ripples are seen in the Cen- 
taurus cluster [51] and M87/Virgo [24]. Concentric ripple- like 
features are also seen surrounding the center of Abell 2052, 
but current analysis shows that these are consistent with weak 
shocks 32 . 

A key difference between shock and sound wave heating 
is that sound wave dissipation depends on transport coeffi- 
cients and shock dissipation does not [31]. The sound wave 
dissipation rate is sensitive to temperature and would be less 
effective in lower temperature systems. This may be offset by 
a higher frequency of outbursts in cooler clusters [31] . One of 
the important implications of the sound waves is that energy 
can be deposited over 47r steradians, and not only along the 
narrow angle of the radio jets. 



Case studies: Abell 2052 and Abell 262 

Abell 2052 is one of the few clusters that showed evidence 
for radio bubbles using data from an X-ray telescope prior to 
the Chandra / XMM-Newton era {ROSAT [52]). It was first 
observed with Chandra for 37 ksec in Cycle 1 in 2000, and 
subsequently observed for 129 ksec in Cycle 6 in 2006 and 497 
ksec in Cycle 10 in 2009. A three-color image from the total 
657 ksec cleaned dataset is shown in Fig. 3, where the col- 
ors represent different energy ranges (red 0.3 — 1.0 keV, green 
1.0-2.0 keV, and blue 2.0-10.0 keV). Clear bubbles are seen 
to the north and south of the cluster center, and a filament ex- 
tends into the northern bubble. Surrounding the bubbles are 
bright shells of emission, and outside of these shells a jump is 
the surface brightness is visible extending around the cluster 
center. A second jump is seen outside this one, and it appears 
sharper to the NE than to the SW. It is possible that this 
feature is a cold front rather than a shock. 

An unsharp-masked image for the total 657 ksec dataset 
is shown in Fig. 4 in grayscale with 1.4 GHz radio emission 
from the Very Large Array Faint Images of the Radio Sky 
at Twenty-cm (VLA FIRST) survey [53] superposed in green 
contours. Concentric ripple-like features are seen surround- 
ing the cluster center, and these are modeled as shocks in [31] 
based on the earlier 163 ksec dataset. These features are at 31 
and 46 kpc from the AGN and the density jumps associated 
with them are both consistent with shocks with Mach number 
1.2. The separation of the potential shock features is consis- 
tent with a radio source outburst rate of 2 x 10^ yr. This 
same cycle time is found when considering the locations of 
outer cavities in the system. A small cavity appears bounded 
by a narrow filament to the NW of the N bubble, and the S 
cavity appears split into two cavities, with the outer one to 
the SE. Assuming the outer cavities rose buoyantly to their 
current positions, we find a cycle time of 2 — 4 x lO'^ yr de- 
pending on whether they rose at 0.5 times or 1 time the sound 
speed. 
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In [32j . we found a mass deposition rate of M = 55±4M0 
yr~^ and a cooling luminosity of Lcooi = 5.4x 10*'^ erg s~^. We 
find the rate of energy input from the buoyantly rising bubbles 
assuming they are filled with relativistic plasma (7 = 4/3) to 
be 3 — 6 X 10*^ erg s~^ depending on the buoyancy rise time. 
The input from shock heating is only 1 x 10'''^ erg s~^, how- 
ever, the combination of heating from shocks and buoyantly 
rising bubbles can offset the cooling in this system. 

Abell 262 was first observed by Chandra in Cycle 2 in 2001 
for 30 ksec. It was subsequently observed in Cycle 8 in 2006 
for 112 ksec. A mass deposition rate of M = lOlgM© yr~^ 
and a cooling luminosity of Lcooi ~ 1-3 x lO''"^ erg s~^ were 
found in [21] . Based on the one clear bubble seen in the earlier 
dataset, it was found that the radio source energy input fell 
more than an order of magnitude short of offsetting the cool- 
ing [21], and it was concluded that the AGN may have had a 
more powerful outburst in the past so that cooling could be 
offset on average. In this study, a repetition rate of 1 x 10* yr 
was assumed for the AGN outbursts. 



Analyzing the longer dataset, however, we found in [22] 
that the radio source can offset the cooling. This analysis was 
based on both a deep Chandra observation and low-frequency 
radio observations. Several more bubbles were found, and a 
tunnel made up of multiple bubbles rising with short time 
intervals was found to the west of the cluster center. Using 
the multiple sets of bubbles, the radio source comes within a 
factor of two of quenching the cooling. Chandra X-ray resid- 
ual images superposed with radio emission at 1440 and 610 
MHz are shown in Fig. 5. It may well be that further analysis 
of deeper multi-frequency observations of objects in samples 
such as [45] that seem to have insufficient energy input rates 
would reveal that the AGN are indeed powerful enough to 
offset cooling. 
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Fig. 1. Composite X-ray (blue, Chandra 657 ksec image), optica! (red. Digitized Sky 
Survey), radio (green/yellow, 1.4 GHz [53]) 5'. 6 X 4'. 7 (230 X 190 kpc) image of Abell 
2052. 
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Fig. 2. Composite X-ray (blue, Chandra 139 ksec image), optical (red, Digitized Sky 
Survey), radio (green/yellow, 1.4 GHz [22]) 5'. 6 X 4'. 7 (110 X 90 kpc) image of Abell 262. 
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Fig. 3. Three-color Chandra 657 ksec 4'. 3x3'. 7 (180 X 150 kpc) image of Abell 2052 
(red = 0.3-1 keV, green = 1-2 keV, blue = 2-10 keV). 




Fig. 4. Unsharp-masked 657 ksec Chandra 2' .S X 2' A (115 X 95 kpc) image of Abell 
2052 with 1.4 GHz radio contours from the VLA FIRST survey superposed. Multiple sets of 
bubbles are seen, as well as ripple-like jumps in surface brightness that are consistent with weak 
shocks. 
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